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ABSTRACT 

We report on a Chandra line spectrum observation of the supergiant, £ Orionis (09.7 lb). 
A 73.4 ks HETGS observation shows a wide range of ionization stages and line strengths over the 
wavelength range of 5 to 26 A. The observed emission lines indicate a range in temperature of 2 
to 10 MK which is consistent with earlier X-ray observations of C Ori. Many lines are spectrally 
resolved showing Doppler broadening of 900 ± 200 kms -1 . The observed He-like ions (O VII, Ne 
IX, Mg XI, and Si XIII) provide information about the spatial distribution of the X-ray emission. 
Although the observations support a wind distribution of X-ray sources, we find three conflicting 
results. First, line diagnostics for Si xni indicate that this line emission forms very close to the 
stellar surface, where the density is of order 10 12 cm~ 3 , but the velocity there is too small to 
produce the shock jump required for the observed ionization level. Second, the strong X-ray 
line profiles arc symmetric and do not show any evidence of Doppler blue-shifted line centroids 
which are expected to accompany an outwardly moving source in a high density wind. Third, 
the observed velocity dispersions do not appear to correlate with the associated X-ray source 
radii velocities, contrary to expectations of wind distributed source models. A composite source 
model involving wind shocks and some magnetic confinement of turbulent hot plasma in a highly 
non-symmetric wind, appears to be needed to explain the line diagnostic anomalies. 
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The Orion Belt supergiant £ Ori (09.7 lb) is 
one of the most thoroughly studied hot stars. It 
has been previously observed in X-rays with the 
Einstein, ROSAT, and ASCA satellites. Observa- 
tions with the Einstein Solid State Spectrometer 
(SSS) found it to be the first O-star to show line 
emission at Sixin and Sxv (Cassinelli & Swank 
1983). ASCA observations show that £ Ori may 
have a weak high temperature component greater 
than 10 MK (Waldron, Corcoran, & Kitamoto 
1997), but the majority of the X-ray emission is 
soft, peaking below 1 keV. The observed X-ray 
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emission is consistent with the standard O-star re- 
lation: L x ~ 10~ 7 L Bo i (Berghofer et al. 1997), 
and has remained relatively stable except for a 
short term of ~ 15% variability observed during a 
ROSAT monitoring study by Berghofer & Schmitt 
(1994). 

It is generally believed that the hot star X-ray 
emission arises from shocks that develop in the 
outflowing stellar wind owing to instabilities asso- 
ciated with the line driven wind mechanism (Lucy 
& White 1980; Owocki, Castor, & Rybicki 1988; 
Feldmeier 1995). The original X-ray models for 
hot stars proposed that the X-rays arise from a 
thin coronal structure located at the base of the 
stellar wind (Cassinelli & Olson 1979; Waldron 
1984). However, the base coronal model was con- 
sidered inappropriate because early observations 
did not show significant attenuation of soft X-rays 
at the oxygen K-shell edge which was expected due 
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to the overlying high density cool wind (Cassinclli 
& Swank 1983). 

The X-ray emission lines produced by a stel- 
lar wind distribution of outwardly propagating 
shocks should show broad profiles with a shape 
that depends on the wind optical depth (Mac- 
Farlane et al. 1991). For an optically thin wind, 
the predicted profiles are flat-topped and symmet- 
ric, whereas, for an optically thick wind, the pro- 
file becomes asymmetric, showing a blue-shifted, 
roughly triangular-shape, owing to the larger at- 
tenuation of the line X-rays from the far side of 
the stellar envelope. The net degree of asymme- 
try is determined by the number of shocks and the 
column density, N w , of the overlying stellar wind. 
The blue- ward extent of the line reflects the largest 
outward velocities of the shock fronts. 

The Chandra spectrometers have the necessary 
spectral resolution and sensitivity to test these 
shock model predictions. In this Letter we present 
a preliminary analysis of the spectral lines from 
our Chandra High-Energy Transmission Grating 
Spectrometer (HETGS) observation of ( Ori. 

2. Observations and Analysis 

Our Chandra HETGS observation of £ Ori was 
performed in two parts: on 2000 April 8 (04^08™ 
- 21 h 37 m UT), and on 2000 April 9 (20 /l 38 m - 
0l h 05 m UT), with effective exposure times of 59.7 
ks and 13.7 ks respectively. £ Ori is the primary 
star of a widely separated binary system with an 
early B-star giant companion, and with Chandra 
the two components are cleanly separated. Here 
we use the standard pipeline processed data from 
our combined 73.4 ks observation to discuss the 
X-ray line spectrum of £ Ori A. The high-energy 
grating (HEG) and the medium-energy grating 
(MEG) spectra are shown in Figure 1. The range 
of ions identified indicate an X-ray temperature 
range of <~ 2 to 10 MK. The triads of He-like lines, 
commonly known as the 11 fir lines" (see Sec 2.1), 
of O VII, Ne ix Mgxi and Sixm are identified in 
Figure 1. 

In this Letter, we focus on the two most ob- 
vious features observed in the HETGS spectrum 
of £ Ori; (1) the unusual line emission properties 
of He-like ions, and; (2) the unexpected combina- 
tion of large line widths and negligible line centroid 
shifts. 
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Fig. 1. — Co-added I s * order, background sub- 
tracted, Chandra HETGS MEG (top panels) and 
HEG (bottom panel) count spectra of £ Ori. The 
ions responsible for the strongest line emissions are 
identified. The bin size is 0.01 A. 

2.1. Analysis of He-like Spectral Line 
Transitions 

Gabriel & Jordan (1969) were the first to 
demonstrate that the fir (forbidden, intercom- 
bination, and resonance) emission lines of He-like 
ions can provide very useful diagnostics for the 
study of X-ray emitting plasmas. They found 
that the line flux ratio 1Z = f/i is density sensi- 
tive, and the ratio Q = (/ + i)/r is temperature 
sensitive. Although these fir line diagnostics have 
been used extensively in solar studies to derive X- 
ray electron densities (n e ) and temperatures (T e ) 
(e.g., Wolfson et al. 1983), the application to O- 
star spectra requires special care. Blumenthal, 
Drake, & Tucker (1972), showed that the presence 
of an external strong UV radiation field can lead 
to overestimates of the density, since not only col- 
lisions, but radiative excitations as well can lead 
to a decrease in the strength of the / line relative 
to the i line. Hence, for O-stars, the observed 1Z 
values for several He-like ions, may no longer lead 
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to a direct measure of n e , but instead provides 
information on the radial distance of the X-ray 
source from the star. 

Using the formalism developed by Blumcnthal 
ct al. (1972), we calculate the radial dependence 
of 1Z in the envelope of £ Ori using a photosphcric 
UV flux model at T eff = 30,000K (Kurucz 1993). 
The radial dependent wind density is derived from 
the stellar and wind parameters of Lamers & Lei- 
thercr (1993) using a standard wind velocity law 
that is modified below 1.02 stellar radii to pro- 
vide a smooth density transition to the photo- 
spheric structure. The radial dependence of the 
UV mean intensity J v {r) is assumed to be geo- 
metric dilution of the stellar radiation, with no 
attenuation, i.e., J v {r) = 4 W(r)H v (T e ff), where 

W(r) = 5 (l - V 1 - ( R *A") 2 )> and H» is S iven 
by Kurucz (1993). The predicted 1Z dependencies 
for Ovil, Ne ix, Mgxi, and Sixin are shown in 
Figure 2. 
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Fig. 2.— The Ovil, Ne ix, Mgxi, and Sixm f/i 
ratio dependence on geometric dilution (radius) of 
the UV radiation field and electron density (ver- 
tical lines) for the circumstellar environment of 
£ Ori. The solid lines include the effects of a UV 
radiation field (appropriate for hot stars), whereas 
the dashed lines neglect these effects. The range 
in the observed f/i ratios are indicated by heavy 
solid lines. The heavy dashed line for Mgxi is an 
extension owing to possible overestimates of the 
/ line strength because of contamination from a 
Mgx line. The corresponding wind densities asso- 
ciated with the radius scale are shown at the top 
of the panel. 



For comparison, we also show the 1Z depen- 
dence for the collision only case, appropriate for 
very low UV flux. 

We find that the predicted 1Z ratios for the 
Chandra observed He-like ions, except Sixm, are 
controlled by the strength of the UV radiation 
field, independent of density, as indicated in Fig- 
ure 2 by the large differences between the colli- 
sional and radiative curves. However, for Sixm 
the predicted 1Z ratio shows three regions of dom- 
inance, collisional (n e > 10 13 cm -3 ), constant 
geometric dilution (< 1.1R*), and a decreasing 
geometric dilution for large radii. This differ- 
ent 1Z dependence for Sixm, as compared to the 
other three He-like ions, is related to the differ- 
ent strengths of the UV radiation field on the two 
sides of the Lyman jump. In the case of Sixm the 
photo-excitation between the upper levels of the / 
and i lines (2 3 S — ► 2 3 P) is produced by radiation 
shortward of 912 A, where the emergent stellar 
flux is <~ 8 times smaller than that longward of 
912 A. Fi gure 2 shows that Sixm marginally al- 
lows for a direct estimate of n e . Since the other 
He-like ions have their excitation transition long- 
ward of 912 A, the observed 1Z ratio can only lead 
to a measure of the geometrical dilution factor (or 
radius). 

Figure 2 shows the measured range of 1Z val- 
ues (given as the heavy solid lines) as a func- 
tion of radius. The 1Z values for Ovil and Ne ix 
are extracted from the MEG spectrum, whereas, 
for Mgxi and Sixm, the HEG spectrum is used. 
From Figure 2 we can read off the radii of the X- 
ray emission line regions, and the associated model 
dependent electron densities. The most interest- 
ing result is the observed 1Z ratio for Si xni of 
1.2 ± 0.5 which indicates that this line emission 
originates from the base of the wind, where the 
electron density is of order 10 12 to 10 13 cm -3 . The 
Mg xi 1Z ratio gives a range in the formation radius 
of 3 to 5 R*, and n e of order 10 9 cm~ 3 . The heavy 
dashed line for Mgxi in Figure 2 represents the 
possibility that the / line may be contaminated 
by a Mgx line since the HEG Mgxi / line shows 
evidence of a red-shift of ~ 0.015 A. The other two 
ions, Ovil and Ne ix, have very weak / emission 
lines, hence, their TZ ratios can only provide upper 
limits on the source radii of 12 and 6 R*, cor- 
responding to lower limits on n e of 4 x 10 7 cm~ 3 , 
and 2 x 10 8 cm -3 respectively. The X-ray location 
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from Mg xi and upper limits from O vn and Ne ix 
shown in Figure 2 suggest that some of the X-rays 
are produced high in the wind, consistent with 
stellar wind shock predictions (Fcldmcicr, Puis, & 
Pauldrach 1997). Furthermore, shocks may still 
be generated above these upper limits, but they 
cannot be the dominant source of the He-like X- 
ray emission lines, otherwise their / line fluxes and 
1Z ratios would be larger than observed. 

The X-ray source radii predicted from the f/i 
analysis (Fig. 2) suggest a systematic ordering of 
the He- like ions, i.e., high energy ions at low radii 
and low energy ions at large radii. In a shock 
picture, a reverse ordering would have been ex- 
pected, since the shock jump temperature should 
be proportional to the local wind velocity. This 
particular anomaly might be explainable by con- 
sidering the radial distance from which the line 
emission can be seen by an external observer. Us- 
ing the wind absorption cross sections of Waldron 
et al. (1998), we find that the optical depth in 
the continuum adjacent to a line reaches unity at 
radii of 1.03, 1.5, 2.8, and 2.2 for Sixm, MgXl, 
Ne ix, and Ovu, respectively. These radii are 
all just less than or equal to the radii predicted 
by our f/i ratio analysis. We suggest, there- 
fore, that our Chandra observations imply that 
the X-rays emission lines originate primarily from 
just above their associated monochromatic X-ray 
" photosphere" where the continuum radial optical 
depth is unity. 

Finally, using the results of Pradhan & Shull 
(1981), Wolfson et al. (1983), and Kccnan et al. 
(1994), we extract T e from the measured Q ratios. 
For Ovu and Ne ix, we find upper limits of 2 and 
6 MK, and for Mgxi and Sixm, ranges of 2 - 11 
and 1-9 MK are indicated. These temperatures 
are consistent with the overall observed ion distri- 
bution illustrated in Figure 1. 

2.2. Analysis of the Emission Line Profiles 

After correcting for instrumental broadening, 
all X-ray lines in £ Ori are seen to be very broad, 
with velocity dispersions in the range 900 ± 200 
kms -1 . Our HETGS observations of ( Ori show 
no evidence of any large Dopplcr blue-shifted 
asymmetric X-ray line profiles. We also find no 
indication that the Gaussian width correlates with 
the ion properties, which suggests that the lines 
form in the same region or at least in regions 



with similar velocity spreads. From the strongest 
HEG lines, we find a distribution in shifts from 
line center of ±0.005 A which is consistent with 
the wavelength uncertainty of the HETGS/HEG 
(Marshall 2000). There are a few MEG spectral 
lines showing blue-shifts of ~ 0.01 A (e.g., Ovill, 
18.97 A), but we also see similar red-shifts. In 
this preliminary analysis we find that there is no 
observed preference for a systematic blue or red- 
shift. 

Our f/i analysis in Sec. 2.1 indicates a distri- 
bution of X-ray sources from <~ 1 to 12 R*. How- 
ever, our observed line profiles appear to show very 
similar line-width and line-symmetry properties, 
which seems inconsistent with their being formed 
over such a large range in radius. To better un- 
derstand this problem, we explore two distributed 
shock models. We only focus on those X-ray lines 
that form above ~ 1.5 R, which arise from shocks 
in the wind (i.e., Mgxi, Ne ix, and Ovu). 

Expanding on the single shock model of Mac- 
Farlane et al. (1991), we developed a shock model 
consisting of 10 spherically symmetric shocks 
equally distributed between 0.4 to 0.97 v^, with 
temperatures ranging from 2 to 10 MK. The shock 
strength scales with the local wind emission mea- 
sure, and the wind absorption cross sections are 
from Waldron et al. (1998). This choice of pa- 
rameters generates an X-ray source distribution 
with radius that is consistent with our f/i analy- 
sis. The model predictions are shown in Figure 3a 
for the two limiting cases, optically thick (maxi- 
mum N w ), and optically thin (N w = 0). In the 
optically thick case (appropriate for a high den- 
sity hot star wind), the wind attenuation reduces 
the contribution to the line emission from the far 
side of the star and the resultant line emission 
tends to arise from the highest density regions 
that are blue-shifted and moving towards us. For 
the optically thin case, the line emission is narrow 
and almost symmetric (with a slight asymmetry 
due to stellar occultation effects) . The narrowness 
arises because emission is proportional to n^, thus 
the strongest line emission occurs from regions 
of low velocity. Comparison of the shock model 
predictions (Fig. 3a), folded through the MEG 
instrumental response, with a typical single ( Ori 
X-ray line profile (Ne x, 12.13 A) is shown in Fig- 
ure 3b, along with the best-fit Gaussian model 
(AA D = 0.042 ±0.003 A). Clearly the optical thick 
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line model does not explain the observed line pro- 
file, and the optically thin case provides a good fit 
to the observed symmetry and broadness. How- 
ever, the optically thin model requires extreme 
conditions. 
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Fig. 3. — X-ray line profile dependence on wave- 
length shift (AA) from line center showing 3 in- 
put line profiles and their expected HETGS/MEG 
count spectra, (a) Model line profiles for a dis- 
tribution of shocks for an optically thick (dashed 
line) and thin wind (dot-dashed line). The best- 
fit Gaussian model profile (solid line) is shown for 
comparison which required a physical velocity dis- 
persion of 1040 ± 75 kms -1 . (b) Comparison of 
the observed Ne X line at 12.13 A (1 sigma error 
bars) with model line profiles folded through the 
instrumental response of the HETGS/MEG. 



The wind could be made thin by increasing 
the ionization structure, which reduces the opac- 
ity (Waldron 1984; MacFarlane, Cohen, & Wang 
1994). However, for normal hot stars, the winds 
cannot be uniformly highly ionized because these 
stars show strong P-Cygni UV line profiles from 
low ion stages such as C iv, which are distributed 
throughout the wind (Lamcrs & Morton 1976). 
An optically thin wind can also be produced by 
reducing the mass loss rate, M since N w is propor- 
tional to M. For example, to maintain an X-ray 



line symmetry comparable to the optically thin 
case shown in Figure 3a, while constraining the 
sources to lie within 12 R», M must be reduced 
by at least an order of magnitude. However, such 
a reduction in M would be inconsistent with wind 
line profile analyses ( Lamers & Morton 1976) and 
free-free continuum observations (Lamers, Waters, 
& Wesselius 1984, Abbott et al. 1980). There 
are other alternatives. For example, a wind with 
a clumpy or non-symmetric structure, has been 
shown to effectively produce an optically thin wind 
even in stars with relatively large M (Moffatt & 
Robert 1994). In addition, if the shocks are pro- 
duced at very large stellar radii where the corre- 
sponding column densities between equal red and 
blue-shifted velocities components are negligible, a 
symmetric line profile can be produced (Feldmeier, 
private communication). However, to achieve a 
difference of no more than 10% between the red 
and blue-shifted velocity components for the Ne X 
line shown in Figure 3b, we find that the shock 
would have to be located at ~ 25 R* which is ~ 
2 times larger than the largest formation radius 
inferred from our f/i analysis. 

3. Discussion 

Our HETGS observation of £ Ori has presented 
several unexpected results. Although the analysis 
of the He-like ions suggest a stellar wind distri- 
bution of X-rays, consistent with a shock model, 
we also find convincing evidence for X-ray emis- 
sion from near the stellar surface. Such a source 
location presents a problem because the wind ve- 
locity would be too small to produce a shock jump 
capable of producing the Si xin ion stage. The ob- 
served X-ray line profiles are very broad and sym- 
metric, and they show no evidence for Dopplcr 
blue-shifted line centroids. Since the observations 
indicate that the shocks must be located within 
12 R*, the observed line profiles imply that a sig- 
nificant reduction in N w or M is required to ex- 
plain the observed line symmetry. Even if we can 
explain the line symmetry problem by applying 
non-symmetric or clumpy wind structure argu- 
ments, we still need an alternative model to ex- 
plain the origin of the high density, near-surface 
X-ray plasma from ( Ori. We suggest that there 
is magnetic loop confinement as was originally pro- 
posed by Cassinelli & Swank (1983) to explain 
their SSS observations of ( Ori. 
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The observed £ Ori surface X-ray plasma den- 
sity is about 2 orders of magnitude larger than 
normal for a quiescent solar magnetic loop (e.g., 
Vaiana & Rosner 1978), but is comparable to so- 
lar flare densities (Phillips ct al. 1996). A sur- 
face distribution of magnetic loop structures could 
provide a mechanism for producing a collection of 
symmetric blue and red-shifted bulk plasma mo- 
tion, and since both components are near the sur- 
face of the star, their emitted X-rays must pass 
through similar amounts of absorbing material, 
producing symmetric line profiles like those ob- 
served in £ Ori. Although the observed velocity 
dispersions are high, they are comparable to pre- 
dicted solar flare velocity dispersions of order 1000 
kms -1 (Antonucci, Bcnna, & Somov 1996). If 
the high density plasma of £ Ori is residing in a 
collection of surface magnetic loops, we can esti- 
mate the required magnetic field by assuming a 
balance between the magnetic and gas pressures 
(i.e., the plasma (3 is unity). Using an n e of 10 12 
cm~ 3 and a T e of 5 MK leads to a magnetic field 
strength > 180 Gauss. This is comparable to the 
field of ~ 200 Gauss detected on the Bl giant (3 
Cephei (Hcnrichs ct al. 2000). The likelihood 
of such fields on hot stars has been investigated 
by Cassinelli & MacGregor (1999) who found that 
fields in hot main sequence stars generated by dy- 
namo action at the interface between the radiative 
core and convective envelope could rise to the sur- 
face. Thus a magnetically confined plasma on the 
surface of a hot star is certainly plausible. 

Because of the r? e dependence of recombination 
that leads to the He-like ions, the fir diagnos- 
tic of density or of source location tend to yield 
values for the densest region producing the lines. 
This was discussed for the density dominated case 
by Brown ct al. (1991). We find evidence that 
X-rays are being produced both at dense regions 
very near the star, and from regions distributed 
throughout the wind. However, the observed X- 
ray line profiles seem incompatible with their be- 
ing distributed as shocks throughout the optically 
thick wind. A more fundamental question that 
we raise is - why are the velocity dispersions es- 
sentially the same regardless of the X-ray source 
location? Our results provide both new facts and 
new challenges for future studies of hot star X-ray 
emission. 
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